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Integrable Cosmological Models Prom Higher Dimensional Einstein Equations 
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We consider the cosmological models for the higher dimensional spacetime which includes the 
curvatures of our space as well as the curvatures of the internal space. We find that the condition for 
the integrability of the cosmological equations is that the total space-time dimensions are D — 10 or 
D = 11 which is exactly the conditions for superstrings or M-theory. We obtain analytic solutions 
with generic initial conditions in the four dimensional Einstein frame and study the accelerating 
universe when both our space and the internal space have negative curvatures. 
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I. INTRODUCTION 

There has been much attention to the understanding 
cosmology from the superstring or M-theory. The recent 
observation of Type la supernovae and CMB measure- 
ment by WMAP indicates our universe is an accelerating 
universe. In String/M- Theory, there is the no-go theo- 
rem [l[ which is an obstacle to realize the accelerating 
universe. The no-go theorem indicates that the warped 
compactification with the static internal space does not 
give rise to the four-dimensional de Sitter spacetime un- 
der some assumptions. This implies that the warped 
compactification can not lead the accelerating universe 
with the static internal space. 

One way to avoid the no-go theorem is to employ 
the time-dependent internal space Q. pace. If the 
curvature of the internal space is negative in the four- 
dimensional Einstein frame, it has been shown that the 
effective Lagrangian has positive potentials [l]-[7[. This 
positive potential gives rise to the acceleration of the 
four-dimensional spacetime in the four-dimensional Ein- 
stein frame. In the four-dimensional Einstein frame, the 
scale factor of the internal space appears as a scalar field 
with potentials which come from the curvature of the 
internal s S-brane solutions also lead to accelerating so- 
lutions [l]-[lB|. The flux field has a role of a positive 
potential and contributes to the acceleration of the four- 
dimensional spacetime. Scalar perturbations of scale fac- 
tors have been shown in 0] that the eternally accelerating 
universe is realized in the eleven-dimensional spacetime 
by the external and internal spaces which possess a nega- 
tive curvature. Recently, it has been also shown by scalar 
perturbations of scale factors that the eternally accelerat- 
ing universe occurs in the ten-dimensional spacetime || . 
The fixed point analysis [3(3] also showed that the eternal 
acceleration is realized if two spaces possess a negative 
curvature. 

In general, it is difficult to solve the Einstein equations 
exactly because of the non-linearity of the Einstein equa- 
tions. Therefore, most of the analytic solutions are spe- 
cial solutions with particular initial conditions. However, 
it is desirable to find integrable models for the analysis 
of the initial conditions for our universe. Up to now, few 



integrable models have been found. In p-brane and cos- 
mological solutions which were inspired with String/M- 
theory, it is known that there exist a few classes of models 
whose solutions can be reduced to the Liouville or Toda 
type [ill- [22j]. The Toda equation is integrable and pro- 
vides exact solutions for us. In [16]- 22], the metric has 
two spatial parts whose curvature is flat for all spaces or 
for one of two spaces. However, if both spaces have cur- 
vatures, it is very difficult to solve even vacuum Einstein 
equations exactly. When our universe starts from the 
quantum era, it is natural to expect that we have spatial 
curvatures whose values are comparable with the curva- 
ture to the direction of time. The curvatures with respect 
to the spatial directions can be regarded as potential en- 
ergies whereas the time dependence can be regarded as 
the kinetic energy of the effective actions. If the universe 
started from quantum fluctuations, it is natural to expect 
that the order of the potential energies are the same as 
the order of the kinetic energy. Therefore, it is highly de- 
sirable to find integrable models with spatial curvatures 
both for our space and the internal space. 

The accelerating universe occurs when two spatial 
parts have the negative curvature as shown in 4]-j)]. 
This analysis performed by the perturbation of scale fac- 
tors. In [3Cj the fixed point analysis was performed and 
also showed that the eternally accelerating universe oc- 
curs with two spatial parts whose curvature is negative. 
In this paper, we will try to solve Z?-dimensional vacuum 
Einstein equations with two homogeneous spaces exactly. 

In order to find integrable cosmological models, we will 
adopt a method for solving Einstein equations analyzed 
in [23j| . In the method used in [23[, the Einstein equa- 
tions can be reduced to an analytically mechanical prob- 
lem with one gauge degree of freedom. The gauge degree 
of freedom originates from the choice of the time vari- 
able. By using time variables which include scalar fields, 
the condition of the integrability has been classified. We 
will review this method in the next section and find that 
we are able to obtain general solutions for those vacuum 
Einstein equations. Actually, it will turn out that the 
classification of [23j was not complete and we will show 
that a new type of integrable models is useful for our 
analysis. 
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We will start with the generic space ansatz with total 
dimensions D and our spatial dimensions d. Strangely 
enough, we will find that the condition for the integrabil- 
ity is that the total dimension is ten or eleven, which is 
exactly correspond to the consistency conditions of su- 
perstrings or M-thcory. If the integrable condition is 
satisfied there is an additional conserved quantity aside 
from the Hamiltonian constraint. Therefore the system 
has two conserved quantities for two dynamical variables 
and then the system reduces to the integrable case. 

In the relation between the integrable con- 

dition and the conserved quantity was investigated from 
the Hamiltonian viewpoint. The same integrability con- 
dition was obtained and it was shown that there exist the 
conserved quantity under the integrability condition. 

The integrable system does not necessarily have the 
simple analytic solutions [13]- [29|. In this paper, it will 
be shown that we are able to derive the analytic solutions 
by the particular choices of the time variable in D = 10 
and D = 11. It is very important to note that the time 
variable used in this paper easily realizes the analytic 
solutions. The other choices of the time variable make it 
difficult to solve the equations of motion. 

We investigated the cosmological behavior in D = 10 
because this case includes the four-dimensional spacetime 
and the six-dimensional internal space. It was found that 
the accelerating universe occurs if two spaces have a neg- 
ative curvature. 

This paper is organized as follows. In the next section, 
we will construct the effective action which has two po- 
tential terms arising from curvatures of the two spaces. 
We will also show that the system is integrable when 
the total dimensions of the space-time is ten or eleven. 
In section three, we will solve Einstein equations com- 
pletely and general solutions will be given. In section 
four, we will consider the cosmological property of the 
case of (D = 10, d = 3) in detail and the analytic so- 
lutions representing the accelerating universe when two 
spatial parts have negative curvature. The asymptotic 
behavior of the metric is also discussed. Finally, we con- 
clude with a discussion of our results and some implica- 
tions. 



dimensional spacetime is given by, 
ds 2 D = G M N(X)dX M dX N 



where 



-e 2n ^dt 2 + e 2A ^ 9ij (x) dx l dx 3 

+ e 2a Wg AB (y)dy A dy B , (1) 



i, 3 = 1, 2, ••• , d 

A, B = d+ 1, d + 2, • • ■ , D - 1 

M, N = 0, 1, • • • , D - 1. 



e 2A(t) an j e 2a(t) are sca i e f ac tors of the d- and (D—d—1)- 
dimensional spaces whose metrics are gij{x) dx 1 dx 3 and 
g ab{u) dy A dy B . respectively. These scale factors depend 
on time variable t. 

We assume that both two spatial manifolds are homo- 
geneous spaces (Einstein spaces); 

^R ij (g(x)) = (d-l)k (d)9ij (x), (2) 
( D - d -VR A B(9(y)) = (D — d — 2) k{ D _ d _ 1} g AB {y), (3) 

where k( d ) and k', D _ d _ 1 s represent the curvature of Ein- 
stein spaces. For our physical d-dimensional space, we 
assume homogeneous and isotropic space. However we do 
not need the explicit representation of gij (x) and gAB (y) 
to derive the effective action. 

The (d + l)-dimensional Einstein frame is realized by 
the following conformal transformations; 



A(t) — A(t) 
n{t) > n{t) - 



1 



d- 1 
1 

d- 1 



(D 



-l)a(i), 
l)a(f). 



(4) 



In fact, the .D-dimensional Einstein— Hilbert Lagrangian 
can be written as 



^GR= y/-g d+1 (x) g D -d-i(y) e dA+r ^ D - d -^ <* R 
^-g d+1 {x)g D - d - 1 (y)e dA+n ( ^R + •••)• 



II. EFFECTIVE ACTION AND CONDITIONS 
FOR INTEGRABILITY 



Under the conformal transformation |4]), we obtain an 
effective action in the (d+ l)-dimensional Einstein frame 

1; 



Let us consider the D-dimcnsional spacetime con- 
structed with two spatial parts which have curvatures. 
For simplicity, we will consider the vacuum Einstein 
Equations. We consider that the D-dimensional space 
consists of two homogeneous spaces whose sizes de- 
pend on time. Namely, the metric ansatz for the D- 



dt 



„dAA 



„dA- 



-A 2 + X 2 a 



i2 } 



C(d)' 



-2A 



k(D-d-i)e d - 



(5) 
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where Vd-i = f dP X X — det (x) yj det gAB (y) and 



X 2 



"(D-d-l) 



(D-d - \){D-2) 
d(d- 1) 2 ' 
(D-d-l)(D-d-2) 
1) 



(6) 



ki 



(D-d-l)- 



The effective action ([5]) shows that the scaie factor of 
the internal space appears as a scalar field with the po- 
tential terms. The two effective potentials are generated 
by the curvature of Einstein spaces which can be easily 
seen from the fact that the potential terms are directly 
proportional to the curvatures. 

We would like to solve equations of motion derived 
from the (d+ l)-dimensional effective action ([5]). In gen- 
eral, Einstein equations are highly non-linear equations 
and difficult to be solved exactly. The effective action 
([5]) indeed results in non-linear equations. A way to an- 
alyze the system is to utilize a gauge degree of freedom 
[2 61 ] . Let us recall that we can choose a time variable via 
the lapse function e n ^> which is a non-dynamical quan- 
tity [26|] . This gauge degree of freedom represents the 
invariance under the coordinate transformation of time. 
Because of this gauge degree of freedom, we can freely 
choose the gauge to solve the system. We will take the 
lapse function as 



e «(*) = e pa(t)+qA(t) 



(7) 



where p and q are any real numbers. We would like to 
look for the more convenient transformation of two dy- 
namical variables. 

We use the following transformation between (A, a) 
and (U+, U-), 



Af 1 + 1 Agg_+1 

e 2A = jj +d -i U J=1 ^ 

(Ml+l) (M 2 + l) 

„2q _ jj(d-l)X jj (d-1) X 



u_ 

d (M]+l) 



(8) 



,2u = u J+ 7=1 jjj 



d (M 2 + l) 
d-1 



where we have defined 
Mi 

M 2 



(d - 2 + q) X + p 
(d-q)X-p ' 
(d - 2 + q) X - p 
(d-q)X+p 



(9) 



By using these variables, we can re-write the effective 
action (0) as [23] 



_ (M x + l)(M a + l) • • n M 1TJ M 2 

i-efi — _ ^ 2 U + U- + A,,/, ( , ( 



+ k(D-d-l) U + 



i I Mj+l 



M 2 + l 



u_ 



{d—l)X 

(10) 



The effective action J dtjC e g(Mx, AI2) can be trans- 
formed to the action of some other parameters Nx,N%, 



J d££ e g(Ni, N2) by the change of time coordinates, dt 
N(£) d£. This is realized by following transformations; 



Ni + l JV 2 +1 



N = V_ 



1 1 "1+1 1 1 "2+1 

L ~<~ M!+l y - L " t "Af 2 + l 



(11) 



These transformations preserve the form of the effective 
action (JTUJ). It is possible to connect a solution in some 
parameters (Mi, M2) to many other solutions by above 
transformations (fTTjl . Because of this gauge degree of 
freedom, we can solve the Einstein equations with a par- 
ticular choice of the parameters. 

In [23j], Mi = Mi = was considered. In our cases, 
this condition leads the following potential 



W = -k 



(d) 



-k. 



-1+7 



U 



If the form of the potential becomes W 



d+ 7 



«{D-d-l) u + 



U+ orW 



Hd) — k(n-d-i) U-, Einstein equa- 
A model with this 



tions are soluble as shown in [2^ 
type of the potential also studied in [24j]. But the above 
potential cannot take such that potential because the 
total dimension has to be D = 1 in order to satisfy 



-1 



d± 



D-1 
(d-l)X 



= 0. Therefore our model is 



not correspond to the model considered in [231 ] - 

A convenient choice of the parameters is Mi = M2 
which implies 
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(d-iy 



■U+U- + k (d) U+U- 



k(D-d-i) U- 



-l+d^T d-^ 



)X 



u_ 



d-1 "^(d-l)X 



(12) 

The above effective Lagrangian shows the second term is 
an interaction similar to a harmonic oscillator and third 
term represents a non-linear interaction which is an ob- 
stacle to solve equations of motion. If the power of U+ 
or U— is simplified, it is possible to solve the equations 
analytically. To perform this procedure, we will impose 
a condition in which the non-linear term in the effective 
action does not depend on U+; 



- 1 



d-1 



d- 



D-2 



= 0, 



(13) 



(d-l)X- 

where X was defined in ©. We will see below that the 
system is intcgrable if this condition is satisfied. Before 
considering the integrability, let us solve the condition 
(fTHjl . The condition (fTB"]) can be rewritten as 



(d- l)(D-d-5) = 4. 



(14) 



From this, we can immediately derive the condition of 
the integrability D and d as follows 



D = 10, d = 3 
£> = 11, d = 2, 5, 



(15) 
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where (0 was used. Note that D = 10 is the critical 
dimension of the superstring theories and D = 1 1 is the 
dimension of the M-theory! Moreover, d — 3 means our 
spacetime is four dimensions. Therefore, we have inte- 
grable cosmological models for a realistic setup. 

We are going to show that the system is integrable if 
the condition (Ti"3"]) is satisfied. By using the condition, we 
can rewrite the effective Lagrangian and the Hamiltonian 
constraint as 



>C n ff — — 



H = 



(d-iy 



■U+U.+k (d) U + U. 



+ k(D-d-i) U 
-^±^U + U_-k {d) U + U- 



2(d+l)/(d-l) 



- k 



{D-d-1) U_ 



2(d+l)/(d-l) 



(16) 



(17) 



where the second equation represents the total energy 
conservation which can be derived by the variation of 
n(t) in the action ([5]). 

We get equations of motion by the variation with re- 
spect to U+ and U— as follows 



(d-ir 

- k 

jd=w & -- k w 



«(<J) U + 

2(d- 



1) 



d-l 
U- = 0. 



u (d+md-i) = 0> (18) 



(19) 



Because of the equation (TT9"]) . we can easily find the fol- 
lowing conserved quantities; 



(d-iy 



1^2 



(20) 



Since the system has two conserved constants (|17p and 
(|2H)l for two dynamical variables, the total system is clas- 
sically integrable. 

In [27|]-[2||j the integrability was discussed from the 
Hamiltonian viewpoint. It is essential idea that they 
looked not only for functions Poisson-commuting with 
the Hamiltonian H, but also for a function F satisfying 
an equation of the form 

{F, H}p. B . = <j>H 

for some unknown function cj). The Hamiltonian con- 
straint H = indicates {F, H}p b. = 4>H = 0, therefore 
the function F becomes a conserved quantity on this 
Hamiltonian constraint. Using this method, same con- 
sequence (fT5|) was obtained in [2jJ-(29|. In our model, 
e satisfies {e, H}p b. — where we used canonical mo- 



menta P. 



+ 



-(4/(<2-l) 2 )[/_ 



the Hamiltonian (fi"7|) and the Poisson bracket 



-(4/(d-l) 2 )[/. 



+ > 



{qi, 12} 



P.B. 



E 

i=+,~ 



dqi_dq2_ _ dqi dq 2 

dUi dPi dPi dU % 



e commutes with the Hamiltonian and then is a con- 
served quantity. This means that the system becomes 
integrable, because the system has two conserved quan- 
tities for two dynamical variables. If the condition (Ti~3")) is 
not satisfied, e dose not commute with the Hamiltonian 
derived from (fT2|) and {e, ifjp.B. 7^ <I>H. In the next sec- 
tion, we are going to show that it is quite easy to derive 
the analytic solutions by the choice of the time variable, 
(Afi = M 2 = 1). 

We can impose the other type of requirement that the 
non-linear term in the effective action (TT2"]) does not de- 
pend on U+. It turns out that interchanging t/+ and 
U- can be achieved by the replacement d — * D — d — 1 
which means that the interchange of the internal and 
the external space. This symmetry should be present 
because we are just considering the evolution of two 
homogeneous spaces. At the level of the effective ac- 
tion this equivalence results from the re-parametrization 
of the time coordinate. For instance, we will take 



-1- 
in H 



Mi + l 
d-l 



d- 



D-2 



(d-l)X 

and impose — 1 - 



= -1 



m 2 +i 

d-l 
2(d-l) 




d-(D-2)/(d-l)X 

d = 6 
d=5,8 



^ (d-l)X J ~ 1 

which gives 
(21) 



where we have used Qj. The effective Lagrangian is iden- 
tical to (frnj) just by interchanging the two spaces. 



III. GENERAL SOLUTIONS 

In the previous section, we have seen that the Einstein 
equations for two homogeneous spaces are integrable if 
the total dimensions are ten or eleven. In this section, 
we will discuss the case D = 10, d — 3 which is most 
relevant for four-dimensional physics. General solutions 
of (D = 11, d = 2) and (D = 11, d = 5) are shown in 
Appendix [A] 

In this case, the equations of motion (fT5T) - (| 1 9 [) and the 
Hamiltonian constraint (|17p are written by 



U + + k (3) U + + A\ 6) U 3 _=0, 

U- + fc (3 ) U- = 0, 

{/+{/_ + fc (3) U+U- + fc (6) ut = 0, 



(22) 
(23) 
(24) 



where fc( 6 ) = 5 k'^ . Let us first consider the equation of 
motion (f2"3")l . The solution of (|23[) can be easily obtained 
as 



~t + A*], 



A\ cos[ y 

r (A 1 t + A 2 , 

Ai cosh[ v / -fc( 3 )t + A 2 ], 



(fc(3) > 0) 

(fc(3) = 0) 
(*<3) < 0) 



(25) 



where A\ and A2 are constants of integrations. These 
equations show that the behavior of U- is controlled by 
the curvature of the three-dimensional Einstein space. 
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Substituting this J7_ into the equation of motion ([22]) . 
we obtain following equations of motion, 



fc (3) > ; U+ + fc (3) U+ 

+ 4k (6) ^A 1 cos[^k^ ) t + A 2 }\ — u . 

, 3 



>(3) 



0; l7+ + 4fc (8) (i4 1 t + A 



(^it+A 2 y 



0, 



fc (3) < ; 17+ + fc (3 ) U+ 



+ 4 k (6) ( A x cosh[ ^-A: (3) t + A 2 ] \ = 0. 



(26) 



These equations of motion concretely show that the U+ 
have received the forced power from the U-. The inter- 
nal space gives the effect of the forced oscillation to the 
motion of U+. 

In the cases of fc( 3 ) < or fc( 3 ) > 0, it is useful to adopt 
the following transformations; 



riven by 
fc (3) > ; E7 4 



B x cos[y'fe(3) t + B 2 
3fc (6) A? 



2,A 



(3) 



t sin[ \l fc( 3 ) t + A 2 



+ cos[ 3(^ (3) t + A 2 )}, 



8 k, 



(3) 



fc( 3) = ; f/ 4 
fc (3 ) < ; E7 4 



Bit + S 2 ~^(A 1 i + A 2 ) b , 



Bi sinh[ ^-fc( 3 ) f + B 2 ] 



fc (6) A? 



i sinh[ \/ —k(3) t + A 



+ J|-icosh[3(y^t + A 2 )]. 



(28) 



fc (3) >0; C/+=/(i) 



3fc (6) Af 



i sin[ w fc( 3 ) t + A% 



8k 



(3) 



cos[ 3(^(3) t + A 2 )] 



fc (3) <0; U+ = f(t) 



3k (6) Al 



fc (6) A? 



£ sinh[ W -fe(3) i + ^2 ] 



cosh[3( A /-A:(3 ) i + A 2 )] 
(27) 



Using these relations, we can simplify the equation of 
motion (|26[) as 



/(*) + *(3) /(*) = 0, (*(3)>0) 
/(*) + (-fe(3))/(*)=0. (fc(3)<0) 

Then, answers for this equation are simply 
f{t) = B x cos[ Jk^ t + B 2 ], (fc (3) > 0) 



f(t) = B x sinh[ ^k^ ) t + B 2 ], (fc (3) < 0) 



where Bi and i? 2 are the constants of integrations. Com- 
bining all these things, we finally obtain the solutions 



The term proportional to k^ indicates the resonance, as 
the frequency of the harmonic and forced oscillation are 
identical. 

The Hamiltonian constraint ([24]) gives a constraint on 
four constants of integrations, 



fc(3) > : 



9fc (6) A? 



C (3) 



fc (3) < 



0; B 1 = 0, 
9fc (6 )A? 



fc (3) B x cos[A 2 - B 2 \ = 0, 



(-fc (3 ))Bi sinhL4 2 - B 2 ] = 0. 

(29) 



We shall consider the metric which is given by 



2ri 2 A tt tt 

e = e = U+U-, 



„2a 




U-' 



(30) 



where we have used flS), ([SJ and M x = M 2 = 1. These 
equations shows that the four-dimensional part is the 
conformal metric, e 2 ™ = e 2 " 4 . Therefore, in the four- 
dimensional Einstein frame, ten-dimensional metric is 



ds 2 w 



,-6ar „2A 



e 2A (-dt 2 + g ij {x)dx l dx : >)' 



„2o 



A j ,B 



9AB(y)dy dy 





-3/2 







f7+t/_ (-eft 2 + 5„ (ie) cfcW) 
1/2 



g AB (y)dy A dy B . (31) 
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For fc( 6 ) ^ 0, metric components are 

(fc(3) >0) 

9k w Ai 



2A _ yK (6) 



«(3) 



cos[A 2 - B 2 ] 



i sin[ ,/y {3) t + A 2 



+ - cos[3(y'fe ( 3 ) t + A 2 )] 



cos[Jk( S )t + A 2 ], 



'9fc(6)^j 

8fc (3) 



C0S[ < + f?2 ] 

COsL4 2 - S 2 ] 



i sin[wfe( 3 ) i + A 2 ] 



+ - cos[3(^fe (3) t + A 2 )] 



(*(3)=0) 



„2A 



„2q 



fc (6) 



(Al* + .4 2 ) 5 +B 2 



(6) 
5A 2 



(Ai< + A 2 ) 5 +B 2 



(Aii + A 2 ), 



Ai t + A 2 



(fc(3) <0) 
2A _ - 9 ^(6)^1 



8(-fc (3) ) 



sinh[ y^fc^ f + B 2 ] 
sinhL4 2 - B 2 ] 



t sinh[ \ — fc( 3 ) t + A 2 



-cosh[ 3(^-fc (3) t + A a )] 



cosh[ 4/ -fe( 3 ) i + A 2 ] 



~9fc(6)^j 

8(-fc (3) ) 



sinh[ ^/-fe( 3 ) t + B-2 



sinh[A 2 - S 2 ] 
+ f sinh[ + A 2 ] 

+ - cosh[3( ^/-/co) i + A 2 ) ] 



cosh[ v /-fc( 3 ) t + A 2 



(32) 



IV. COSMOLOGICAL CHARACTERISTIC FOR 
(D=10, D=3) 

The solutions that we have obtained in the previous 
section include a metric which has realistic dimensions 



D = 10, d = 3. The total dimension of this spacetime 
is ten dimension equal to the critical dimension of su- 
perstrings and the physical space-time has four dimen- 
sions. Therefore, it is very interesting how the universe 
evolve with time. In this section, we shall consider the 
behavior of the ten-dimensional spacetime. We will ana- 
lytically show that the four-dimensional part of the ten- 
dimensional spacetime accelerates eternally, which has 
been analyzed in [f| by the qualitative method and in 
[30l | by the fixed point analysis on the phase space. 



We shall consider metric 
la take oscillatory behavior. 



J2 A 



In fc (3) > 0, 



J1A 



and 

starts from zero and 
end up with zero because i/_ = A\ cos[^/fc( 3 ) t + A\\ 
oscillates between two zeros. On the other hand, e 2a 
diverges when U- and e 2A become zero, and then, the 
case of fc(3) > may not have a stable internal space. 

Similarly, in Av 3 ) = 0, the scale factor of the internal 
space diverges at U- — A\ t + A 2 = and e 2 A takes zero 
when U- = 0. For fc( 3 ) = and fc 6 < 0, the asymptotic 
behavior becomes e 2A —> t 6 and e 2a — > t 2 at t — ► 00 
and the ten-dimensional metric (|31j) has the behavior as 
follows 



dsi 



(-dt 2 



ds\ 



)+* 



2 dsl 



(33) 



This means that the four-dimensional part of the met- 
ric does not depend on t in the ten-dimensional frame. 
If we take ds 2 to the three-dimensional Euclid space, the 
four-dimensional part becomes the Minkowski spacetime 
at t — > 00. The internal space becomes large in this re- 
gion. 

If fc( 3 ) < 0, an interesting phenomenon occurs. The 
three-dimensional space expands with acceleration eter- 
nally. This acceleration is extracted from the negative 
curvature of the internal space, fygj < 0. The curvature 
of the internal space acts like a positive cosmological con- 
stant in four dimensions. We assume the internal space 
is the Einstein space with the negative curvature. This 
case, fc( 3 ) < and k(g) < 0, is equivalent to the situ- 
ation suggested in [5] in which it was shown by scalar 
perturbations of scale factors that the acceleration of the 
three-dimensional s pac e occurs in the four-dimensional 
Einstein frame. In [301 ] the fixed point analysis also in- 
dicated that the eternal acceleration occurs for fc( 3 ) < 

and fc(6) < 0. We can confirm these facts by using the 
analytic solutions. 

Defining a proper time dr = ^/U^ U— dt, the veloc- 
ity and acceleration of e A are given by de A /dr and 
d 2 e A /dr 2 . These velocity and acceleration are shown 



in fig.pQ-0 where we neglect overall factor in 
From these figures, we can extract following facts. The 
three-dimensional space evolves with the very large posi- 
tive velocity and negative acceleration, then the three- 
dimensional space decelerates quickly at a first stage. 
The acceleration of the three-dimensional space turns to 
positive at some time, and then the acceleration decreases 
gradually. The three-dimensional space finally expands 
with a positive velocity and the zero acceleration at the 



7 



the four-dimensional frame, 



2.2 




FIG. 1: The velocity of e A . k 



(3) 



-1, A 2 = 0.5, B 2 = 0. 



d 2 e A 
dr* 0.02 

0.01 

-0.01 

-0.02 
-0.03 




dsi 



t-*l 2 [t{~dt 2 +d S \)} 



t^ds 2 . 



(36) 



The proper time for the four-dimensional frame is given 
by t ~ i 3 / 2 and the three-dimensional space expands 
with e 2 
is d 2 e A /dr 
diverges at r — > 

As a final example, we shall consider 



r 2 / 3 . The acceleration of this scale factor 
2 - — (2/9)r _5/ ' 3 and then the acceleration 



and 



fc(3) < 0. In this case, we can find that the constants of 
integrations satisfy A2 = B2 in l|29p. Using 
and ([30]) . the asymptotic behavior of e 2A 



„2 A 



e~" — ► const, and e~" — » e 
ten-dimensional metric leads to 



2 



at t 



dsi 



(-dt 2 + dsi) + dsj. 



The 



(37) 



This metric shows that the internal space does not de- 
pend on t. The proper time is defined as t ~ e V~ fc < 3 ) * 
in ten dimension. Therefore, the ten-dimensional metric 
is represented as 



FIG. 2: The acceleration of e A . 



fe(3) = 



-1, A 2 = 0.5, B 2 = 0. 



infinite future. In particular, the three-dimensional space 
accelerates forever. This fact coincides with Q and [30||. 

We shall consider the asymptotic behavior of e 2 A and 
e 2a in ([32]>. J7_ = Ai cosh[ v / -fc (3 ) t + A 2 ] is not singular 
for appropriate constants of integrations. A\ and A2, 
which implies e 2a may not diverge when e 2A — > 0. For 



£ < 3 > r and 



ds 2 



C (3) 



\-dt 2 +ds 2 3 + ds 2 6 



(34) 



It is found that this ten-dimensional metric is confor- 



mally equivalent to — dt 2 + ds 2 + ds\ at large t. In four- 
dimensional Einstein frame, this metric (1341) behaves as 



ds 2 10 



(-dt 2 + ds 2 3 )} 



+ e V^) t ds 2 . 



(35) 



The above metric ([33]) shows that the internal space 
also expands. The proper time, for the four-dimensional 
frame, is given by r ~ e 2 V~ k ( 3 > 1 and then, the three- 
dimensional space expands with e 2A ~ r 2 which shows 
that the expansion has the uniform velocity at large r. 
The three-dimensional space has a negative curvature 
fc(3) < in this case. Therefore the four-dimensional 
part can be the Milne universe for the four-dimensional 
frame. The metric (1351) also indicates that the internal 

-k( 3 ) t 



space becomes large with e 



>.v- 



and then it 

is intuitively expected that the curvature of the internal 
space to decrease. In fact, it is found in (O that the 
potential term —k^e~ Sa vanishes at t — > 00. 
E t -f 0, ds 2 w -» t- x l 2 (-dt 2 + ds 2 ) 



t^-^dsg and in 



dsi 



(-dr 2 + t 2 ds 2 3 ) + d 4 



(38) 



whose structure is the product space of the Milne uni- 
verse and a flat six-dimensional space. It is possible to 
transform the Milne universe into the Minkowski space- 
time by coordinate transformations. In this case, the 
above metric (|38p becomes the product spacetime with 
the four-dimensional Minkowski spacetime and the flat 
six-dimensional internal space at t — > 00. 



V. CONCLUSIONS 

We have considered the vacuum Einstein equations 
in the D-dimensional spacetime and obtained integrable 
cosmological models. Those Einstein equations have two 
potential terms arising from the curvature of the d- and 
(D — d — l)-dimensional Einstein spaces. It was thought 
that solving Einstein equations with two curved spaces 
are very difficult. However we have pointed out that 
the total dimension should be D = 10 or D = 11 to 
make those Einstein equations integrable as cosmological 
models. The integrability is guaranteed by the conserved 
quantity which commutes with the Hamiltonian. The 
integrable system does not necessarily have the analytic 
solutions. It is very important to note that the time 
variable used in this paper easily realizes the analytic so- 
lutions. It is interesting that models with superstrings 
or M-theory are more tractable as cosmological models 
than other dimensional models. 

For (D = 10, d = 3), we have obtained the accel- 
erating universe with two spatial parts whose curvature 
is negative. The three-dimensional space expands with 
the acceleration, but the six-dimensional internal space 
also expands. The external space finally approaches the 
expansion whose acceleration tends to zero. It may be 
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difficult for this case to give an account of the realistic 
acceleration at late time. To obtain more realistic mod- 
els, we need to construct a model whose internal space is 
fixed dynamically whereas our space is going to expand 
more drastically. In the context of pure gravity solu- 
tions we have treated in this paper, we cannot get such 
solutions. The flux field, dilaton and the world volume 
actions may have a role for the interesting behavior such 
as fixing the internal spaces [251 ]. 

It would be more interesting to find solutions for more 
general setup. 



Acknowledgments 

We would like to thank N. D. Hari Dass for valuable 
discussions and Eiichi Takasugi and Tetsuyuki Yukawa 
for useful comments. We would like to thank Nobuyoshi 
Ohta for informing us papers [27l|-[30l|. 



APPENDIX A: GENERAL SOLUTIONS IN D = 11 

We can obtain other solutions corresponding to (D = 
11, d — 2) and (D = 11, d = 5) from the equations (f]~6|) 
and fT7]) . 

In (D = 11, d = 2), the equations of motion and the 
Hamiltonian constraint are 



4U 

4iL + k 



^(2) u+ 



6k (8) U 5 _=0, 



(2) 



U- = 0, 



4U+U- + k (2) U+U- + k (s) U e 



«(2) 

and for (D = 11, d= 5), 
1 



U 4 



k( 5 ) U-i 



3k (5) U 2 _=0, 



4 

1 • 



f/_ + fc (5) f/_ = 0, 



U+U- + k {5) U+U- + k (5) Ui = 0. 



(5) 



For (D = 11, d = 2), solutions are given by 



A 1 cos[ t + A 2 



U- = iA 1 t + A 2 



A l cosh[^-J^i + A 2 



(fc(2) > 0) 
(fc(2) = 0) 
(fc(2) < 0) 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 
(A6) 



(A7) 



k(2) > ; 



U+ 



- B\ cos[ - 
15kr 8) Al 



C(2) 



t + B 2 ] 



16, fk 



(2) 



t sin[ ■ 



Ht 



(2) 



t + A 2 



15fe (8) Ag 
64fc (2) 



cos[3( ■ 



*(2) 



t + A 2 )} 



64 k 



(2) 



cos[5( ■ 



C(2) 



t + ^ 2 ); 



e (2) = ; U+ = B 1 1 + B 2 



"(8) 

28Aj 



(A!< + A 2 ) 7 , 



fc (2) < ; C/+ = B x sinh[ t + B 2 ] 



15* (8) A5 , sinh[ V^ i + j4i] 



16y^k 



(2) 



15fe/ 8 -)A? r , \/—k(2) 
—^ C osh{3(^lt + A 2 )} 



+ 



(2) 

cosh[5( 



64 fc 



(2) 



t + A 2 )}, 



(A8) 



and constraints are given by 

fc( 2) > ; ( 2 + fc (2) B i C0S t A 2 - S 2 ] = 0, 
fc (2) =0; Bi=0, 



(A9) 



fc(2) < ; + {-k {2) )Bx sinh[A 2 - B 2 ] = 0. 



For (D = 11, d — 5), solutions are 

Ai cos[2^/fc( 5 ) t + A 2 ], 
L7_ = ^ Ait + A 2 , 



(*(6) > 0) 

(fc(B) = 0) 

Ai cosh[2 A /-*(B}t + A 2 ], (%<0) 



(A10) 



fc (5) >0; U+ = Bx cos[2^k {5) t + B 2 ] 
3fc (B) A? 



2 A. 



(■•5) 



2fc, 



(5) 



cos[2(2«/* (B ) t + Aa)], 



fc( 5) =0; (7+ = B 1 t + B 2 -^(A 1 t + A 2 ) 4 , 



fc (5) <0; C/+ = Bi sinh[2^/-fc (5) t + B 2 
3k {5) Aj 



2 k 



(5) 



2k 



(5) 



cosh[2(2y'-fc (5) i + A 2 )]. 

(All) 



The constraints are 



fc( 5 ) > 

fc( 5 ) = 
fc(5) < 



k {5) Ax B x cos[A 2 - B 2 ) = 0, 

Bi = 0, (A12) 

(-k {5) )A 1 B 1 sinh[A 2 - B 2 ] = 0. 
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For (D = 11, d = 2) and three dimensional Einstein 
frame, the metric is 



{U+U-f (-dt 2 + gij (x) dx'dx 3 ) 

1/3 





-8/3 







U- 



9ab(v) dy A dy B , 

(A13) 



the metric is given by 
ds 2 n = ( j±- I [{U+U-) l ' 2 {-dt 2 + gij (x) dx'dx 1 ) 



U. 



2/3 



jf- I 9AB(y)dy A dy B , 



(A15) 



1/3 



=e^" - {U+U-f, 
For (D = 11, (i = 5) and six dimensional Einstein frame 



h-) • ( A14 ) e 2A = e 2 ^(U + U^ 2 , 



2/3 



(A16) 
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